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I. INTRODUCTION 

The differential cross section of elastic scattering of p+p collisions allows one to study the internal structure 
of protons using the theory of diffraction. Varying the momentum transfer one can change the resolution of the 
investigation: increasing the momentum transfer corresponds to looking more and more deeply inside the structure 
of protons. One of the fundamental outcomes of diffractive p+p scattering studies was the indication that protons 
have a finite size and a complicated internal structure, thus the protons can be considered as composite objects. 

Our interest in this problem has been triggered by two factors: an interesting series of recent theoretical work and 
also new data from the TOTEM experiment at CERN LHC. These are detailed below. 

Recently, we became aware of an inspiring series of papers of Bialas, Bzdak and collaborators, who studied elastic 
proton-proton [l[ , pion- proton Q and nucleus-nucleus collisions [1, H| in a framework where the proton was considered 
as a composite object that contains correlated quark and diquark constituents. In this work, we confirm their main 
conclusion: the quark-diquark model of nucleon structure at low momentum transfer does capture the main features 
of this problem and indeed it deserves a closer, more detailed attention. The first results that form the basis of our 
current study were made public in a proceeding material [j| . In that study we already reported about several details of 
the present investigation. In contrast to the paper of Bialas and Bzdak in ref. [l| , we not only included an estimation 
of the best values of the model parameters, but also determined their errors and also the fit quality. 

In this paper, we refined our first results [B| by decreasing the number of fit parameters from five to three which 
leads to a far better error estimation, and less correlations among the remaining three model parameters. In order 
to reach these goals, we utilized standard experimental techniques, such as multi-parameter optimalization or fitting 
with the help of the MINUIT function minimalization and multi-parameter optimalization package [6J. 

In addition to these simple, straightforward and interesting theoretical investigations of elastic scattering data from 
CERN ISR in the energy range of \/s — 23.5, 30.7, 52.9 and 62.5 GeV, that were already analyzed in ref. [l], new 
clastic scattering data became available recently at y/s = 7 TeV [7] from the CERN LHC experiment TOTEM. The 
new data triggered a wide s pec trum of theoretical investigations, includin g Re gge theory based studies involving 
perturbative QCD(BFKL) [12[, as well as study with hidden dimensions [13| . Extrapolation to the domain of cosmic 
rays [l4| become available. We have tested both variations of the model of Bialas and Bzdak of ref. not only at 
ISR energies but also at the currently available highest LHC energies on recent TOTEM data, in order to learn more 
details about the evolution of the properties of p+p elastic interactions in the recently opened, few TeV energy range. 
The two variations of the Bialas and Bzdak model allowed us to obtain an effective radius parameter i? e ff which is 
model independent within errors. A simple relation between the total cross section and the effective radius is also 
presented in the discussions. 

Most of the arguments for the composite structure of hadrons have been derived from the studies of lepton-hadron 
interactions. The standard picture is that hadrons are either mesons, composed of valence quarks and anti-quarks, 
or (anti)baryons, composed of three valence (anti)quarks, that carry the quantum numbers, while the electrically 
neutral gluons carry color charges and provide the binding among the quarks and anti-quarks. The exact contribution 
of quarks, gluons, and the sea of virtual quark-antiquark pairs and gluons to certain hadronic properties e.g. spin 
is still under detailed investigation. Also, more than 10 exotic hadronic resonances called X, Y and Z states were 
recently discovered in electron-positron collisions at the world's highest luminosities in the BELLE experiment at 
KEK. These hadronic states cannot be interpreted in the standard picture of quark-antiquark or three (anti)quark 
bound states, according to refs. [HI, [l6j]. Thus even nowadays there are still several open questions that are related 
to the compositeness of the hadrons in general. As gluons do not interact directly with leptons, their properties are 
best explored with the help of the strong hadronic interactions. For example, the gluon contribution to the proton 
spin is investigated with the help of polarized proton - polarized proton collisions at RHIC [175, but it is still not 
fully constrained. In this work, we focus on the effects of internal correlations between the quarks inside the protons, 
reporting on a detailed study of proton-proton elastic scattering at several ISR and also at the currently top available 
LHC energies. Let us recall, that a similar analysis involving three independent quarks was not able to properly 
describe the ISR data 18]. In that model quarks were considered as "dressed" valence quarks in the sense that they 
contain the gluonic and qq contribution as well, as if the glue would be concentrated around pointlike valence quarks. 
Thirty-five years after the three independent quark model of ref. [181 ] another three-quark model of the protons was 
proposed in ref. [l|, that included interesting correlations between two dressed valence quarks to form a diquark. In 
the present investigation this quark-diquark model [J . is compared to data in details. 

The quark-diquark picture of elastic p+p scattering is similar to the Glauber optical model [r| in nuclear physics. 
This Glauber model, developed originally for nuclear multiple scattering problems like cross sections of protons 
and neutrons on deuteron, became a standard model of high energy interactions in nuclear physics where multiple 
interactions are built up from superpositions of nucleon-nucleon scattering. This model became a fundamental and 
successful tool in describing nuclear collisions at high energy po| . 

The body of this manuscript is organized as follows: the theoretical model of Bialas and Bzdak is recapitulated and 
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summarized in Section [TH including two cases. In the first case, the proton is modelled as a quark-diquark composite 
system, while in the second case, the diquark is assumed to have a quark-quark internal structure. Section lTITl contains 
our final results, which is based on a multiparameter optimalization procedure, utilizing the CERN MINUIT package 
6]. Finally, we summarize and conclude in Section [IVI 



II. ELASTIC SCATTERING IN THE QUARK-DIQUARK MODEL 



We describe proton-proton interactions as collision of two systems, each one composed of a dressed quark and 
diquark. The p+p collision is schematically illustrated on Fig. [TJ 

^\ Proton 1 




FIG. 1. Scheme of the scattering of two protons, when the proton is assumed to have a quark-diquark structure and the 
diquark is assumed to be scattered as a single entity. This is just a snapshot and all the model parameters follow a Gaussian 
distribution. Note, that a center of mass energy dependent Lorentz-contraction determines the longitudinal scale parameters. 



The interaction between quarks and diquarks is assumed to be purely absorptive. Consequently the amplitude has 
no real part and the imaginary part - dominating at high energy - is given by the absorption of the incoming particle 
wave, represented by the inelastic (non-diffractive) collisions. 

In the impact parameter space the inelastic proton-proton cross-section for a fixed impact parameter b can be given 
by the following formula [1] 

+oo +oo 

<r(6) = / ... / d 2 s 9 d 2 s'd 2 s d d 2 s^(s^,s^)Z?(s^\sdV( s ^ 5 s d;s^',Sd';&), (1) 



where iq, sT q ' and s*d, s*d are the transverse positions of the quarks and diquarks respectively. The integrand is 
a product of quark-diquark distributions of the incoming protons and a er function which gives the probability of 
inelastic interaction at given impact parameter b and at given quark, diquark positions. 

The quark-diquark distribution inside the nucleon is taken into account with the following Gaussian 

D (*9. = l -^rz~ {s2 -> +sl)/Kd 5 2 (sd + Xa- q ), A = m q /m d , (2) 

where R q d is the RMS of the separation between the CMS of the diquark and the remaining quark in the proton. A 
is the mass ratio of the quark and the diquark. Obviously 1/2 < A < 1, where 1/2 would indicate a loosely bound 
diquark. The two dimensional delta function fixes the position of the center of mass in the transverse plane. 

Elastic interactions are independent, accordingly the probability distribution of elastic proton-proton collision is the 
product of the probability distribution of elastic interactions of their constituents (2lL |22| 

cr(s q ,s d ;s q ',s d ';b) = 1 - [l - cr ab (b + s„' - s&') . (3) 

a.b£ {q,d} 
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The inelastic differential cross-sections are parametrized with Gaussian distributions 

Tab (si = A ab e- s2 / R ^, R 2 ab = R 2 a + Rl (4) 

where R ab is the variance of having an inelastic collision, which is calculated from the sum of the squared R q , R d 
radius parameters; the A ab parameters are the amplitudes. From unitarity the elastic amplitude in impact parameter 
space 



t et (b) = 1-^1 -a(b). (5) 

As it was mentioned the real part of the amplitude is ignored. Recently, the important role of the real part of the 
elastic scattering amplitude in shaping da/dt at the dip and in the Orear region was highlighted in (23l. l24j|. 

The elastic amplitude in momentum transfer representation is the Fourier-transform of the amplitude in impact 
parameter space 

+ OO +00 +OG 

T(A) = J J t el (by Krb d 2 b = 2ir J t el (6) Jo (Aft) bdb, (6) 

— oo — oc 

where A = |A|, b = \b\ and Jo is the zeroth Bessel- function of the first kind. Then the elastic differential cross section 
reads as 



A. Model p — (q,d): The diquark is assumed to scatter as a single entity 

The subject of this section is to analyze, as the first investigated case, that variant of the model of Bialas and 
Bzdak, when the quark and the diquark is considered to scatter as one entity as indicated on Fig. Q] In this case, the 
number of free parameters can be reduced if we assume that the number of partons is twice as many in the diquark 
than in the quark. From the inelastic differential cross sections (g]) the total inelastic cross sections are 

a ab = J J Tab (s)d 2 s = TrA ab R 2 abl a,b € {q, d}. (8) 

— oo — oo 

Our assumption tells us that 

<Jqq ■ <Jqd <J dd = 1 \ 2 : A, (9) 

from which we can deduce the following expressions 



4i?2 4i?2 
A q d = A qq ^ 2 — ^ 2 , A dd — A qq -j^-, (10) 

which means that every A ab parameter can be expressed in term of A qq . With these ingredients the calculation of each 
term in (TTJ reduces to Gaussian integrations. Two of the Dirac 5 functions in (fT]) induce the following transformation 
in the transverse diquark and quark position variables 

Sd = -As^, S d ' = -\s q . (11) 

Hence four Gaussian integration remain, which lead us to the following result [l| 

+00 +00 

J J d 2 s? d 2 4 e ~ 2 H^+4 2 ) e "^( b -^+<) 2 e ~^("-^+^) 2 (12) 

— 00 —00 

x e -c dq (b-s d +s' q y e -c dd (b-s d +s' d ) 2 _ ^_ e ~b 2 ^^ 
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where the coefficients c ao are abbrevations, and 



4v + (l + A) {c qd + c dq ) 



\v + c qq + A c dd l 



(13) 



+ (1 - A) V {C q d + C dq ) + (1 + A) CqdCdq 



while 



r = 



4V + (1 + A) 2 (C qd + C d q) V (c qq + C d d) + (1 + A) 2 C qq C d d 



AV + (1 + A) 2 (Cqq + C d d) V {c q d + C d q) + (1 + A) 2 Cq d C d q 



(14) 



B. Model p — (q, (q,q)): The diquark is assumed to scatter as a composite object 



When the diquark is assumed to scatter as a composite object, that includes two valence quarks, the scheme of 
elastic p+p scattering is illustrated on Fig. [5J Following Bialas and Bzdak [l[, the quark distribution inside the 
diquark is supposed to have the following Gaussian shape 

D (Sql,Sq 2 ) - ^e-W^V 2 ^ 2 (Sq X + S q 2 ) , (15) 

where Sq\ and Sq2 are the transverse quark positions inside the diquark, and 

d 2 = B 2 d - R 2 q (16) 
is the RMS of the separation of quarks inside the diquark, calculated from the diquark and quark radius parameters. 





FIG. 2. The scattering situation of the two protons when the diquark is assumed to be composed of two quarks, and the proton 
symbolically can be written as p=(q,(qq)). This is just a snapshot and all the model parameters follow a Gaussian distribution. 



If the internal structure of the diquark is given by eq. (1151) then the a q d, <Jdq and add inelastic differential cross 
sections eq. (|4]) can be calculated from a qq using an expansion analogous to expression eq. (J3)). The result for a qd is 
the following 

- ^T^"^ " 4f e ~ s2 ' K ' (17) 

add is a bit more complicated [![. The inelastic cross section of eq. (fTJ) for this model is again obtained using 
expression (TH21) . 
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In their original paper, Bialas and Bzdak [l| fixed the model parameters by the total cross section, the slope 
parameter B, the position of the dip and the position of the first diffractive maximum after the dip. This way, they 
determined the best values of the model parameters essentially by the method of numerical solution of four equations 
containing four parameters, and found, that the resulting solution gives a good overall description of elastic scattering 
data at ISR energies. 

In the first phase of our studies, we confirmed their results, and noted that if we want to determine the errors on 
the model parameters a different strategy is needed. As we were interested in finding significant dependence of the 
model parameters on the colliding energy, for which the estimation of the errors on the model parameters is essential, 
we started to improve on the original method of Bialas and Bzdak by a multiparameter fit utilizing all the available 
information in the data set, determining the best confidence levels and the errors on the model parameters. 

In the next section, we present the final results of this multi-parameter optimalization method. This manner, we 
improve on the earlier analysis of Bialas and Bzdak |l| with the help of the CERN Minuit package, in both cases: 
when the diquark acts as a single object, and also when the diquark is considered as an object composed of two 
quarks. We perform these analysis for all the energies where p+p elastic scattering data are available from the ISR 
collider, and we apply the same kind of analysis also to recent TOTEM data at s/s = 7 TeV at CERN LHC energy. 

Our final analysis results of p+p elastic scattering data are presented for data sets that are restricted to exactly the 
same kinematic domain at all the five considered energies. We also study, in addition to the differential cross-section of 
elastic scattering, the effect of including the measured total p+p scattering cross-section data to the fitting procedure, 
at all the four ISR and also at y/s = 7 TeV LHC energies. 

Finally, we summarize and conclude. 



7 



III. FINAL MINUIT FIT RESULTS TO P+P ELASTIC SCATTERING DATA AT ISR AND LHC 

ENERGIES 

A. Model p — (q,d): The diquark is assumed to scatter as a single entity 

In this section the MINUIT fit results are presented for the ISR [H, [2(| and TOTEM Q proton-proton elastic 
scattering data considering the scenario when the diquark is assumed to act as a single entity in this scattering 
process. Our preliminary analysis of these data, that did not yet study the correlations between model parameters 
and also did not evaluate what happens if some of the model parameters are fixed, was presented in a recent conference 
contribution Q . 

Our final results are illustrated on Figs. [5][7] The confidence levels, and model parameters together with their 
errors are presented in Table HI The calculated total elastic cross sections, including their uncertainty, were evaluated 
from the MINUIT fits to the differential cross section data, restricted to the same intermediate t elastic scattering. 
At the end of this section we also study and discuss, what happens when the measured total elastic cross sections are 
also added, as an additional data point, to the optimalization procedure. 

The ratios of the inelastic cross sections were fixed with eq. ©, in order to decrease the number of free fit 
parameters, therefore these ratios will be provided only for the case, when the diquark is assumed to be a composite 
object. 

Some additional remarks are due before our final results are presented. The Bialas - Bzdak model shows a singular 
behaviour at the diffractive minimum or dip position, which is due to the lack of a real part in its amplitude. In a 
Minuit fit, such an unphysical fit region may completely dominate the fit results. In order to avoid such an artefact 
and to obtain a meaningful fit result, we have excluded 3 data points from the optimalization procedure, that were 
closest to this dip region. These points are shown in red (color online) or on the plots. We have checked that leaving 
out 5 or 7 points did not change the results. 

Another important remark is that the TOTEM data covers the |t| range from 0.36 GeV up to 2.5 GeV and this 
range is applied in our minimization procedure to allow a comparison between the ISR and TOTEM results. Note that 
Bialas and Bzdak adjusted their model in a different way: they demanded, that four important features of the data 
are described correctly, namely (i) the total inelastic cross-section and (ii) the slope of the differential cross-section as 
extrapolated to t = 0, and (iii) the position of the first diffractive minimum and finally (iv) the height of the diffractive 
maximum after the diffractive minimum [l|. A different strategy is followed here, since we have fitted the theoretical 
curve directly to the experimental data points using the CERN MINUIT package by minimizing x 2 . However, we 
selected the same reasonable /-range for all data sets, and we left out data points around the diffractive minimum, 
where this model clearly breaks down. These data points are clearly indicated with filled (color online red) circles in 
the subsequent Figures, where the fit results are shown and can be compared to data directly. 

In a preliminary analysis, reported in ref. all model parameters were optimized, including also the values of A 
and A qq . Those preliminary results are compatible with the final fits presented here, however, in those results large 
correlations remained among the fit parameters, and their errors were correspondingly large. The optimalization 
process indicated that the parameters A and A qq were in the range of their nominal value and fixing them resulted in 
considerable reduction of correlations among the other model parameters. Hence in the final fits presented here we 
utilized fixed values of A = 0.5 and A ao = 1.0. 



8 



p+p -> p+p, diquark as a single entity at \ s=23.5 GeV 
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FIG. 3. (Color online.) Results of MINUIT fits at ISR energies when the diquark is assumed to scatter as a single entity. Top 
panel shows the data points and the result of the best fit, while the lower panel shows the relative deviation of the model from 
data in units of measured error bars. As the model is singular around the dip, 3 data points that are located closest to this 
difiractive minimum and are indicated with filled (red) circles in this Figure, were left out from the fit. The optimalization was 
restricted to the 0.36 - 2.5 GeV \t\ range, so that a fair comparison could subsequently be made with most recent TOTEM 
results as given in ref. Q • The best fit is indicated with a solid (black) line in this range and its extrapolation to low values of 
t are also shown. The confidence level, after fixing the values of A and A qq , is still higher than 0.1%, which indicates that this 
fit quality is statistically acceptable. 



p+p p+p, diquark as a single entity at Vs=30.7 GeV 



| 10 

B 

E 1 

§10" 1 

o o 

£10 3 



"D 

10 5 



A qq and A. are fixed 
Fit range = 0.36< t < 2.50 GeV 2 
X 2 / NDF = 26.73/32= 0.84 
CL = 73.053 % 



lo- 
ur 7 k 

10' 



10" 

20 
15 
10 

: 5 
o 

-5 
T -10 
jS -15 
S "20 f 




o Data points 

• Not used 
— Fitted theory 
Extrapolation 



nnnnn e e _^ 



— O Q. 



I = 0.50 
A qq = 1.00 

R ^ = 0.34 + 0.01 fm 

qd 

R a = 0.25 +0.01 fm 



R d = 0.71 ± 0.01 fm 



J I I I I I I I L 



J I I I I I I I I I I I I L 



o 



• • 9 • • 
••••• ••••• • • • 



_i i i I i i_ 



_l i i i i I i i i i I i i i 



0.5 



1 2 1-5 2 

-t [GeV 2 ] 



2.5 



FIG. 4. (Color online.) Same as Fig. El but for the energy *fs = 30.7 GeV. 
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p+p -» p+p, diquark as a single entity at \ s=52.8 GeV 
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FIG. 5. (Color online.) Same as Fig. [3] but for the energy y^i = 52.8 GeV. 
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p+p -» p+p, diquark as a single entity at \ s=62.5 GeV 
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FIG. 6. (Color online.) Same as Fig. [3] but for the energy ^/s = 62.5 GeV. 
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p+p -> p+p, diquark as a single entity at \s=7000.0 GeV 
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FIG. 7. (Color online.) The result of the fit at LHC at 7 TeV when the diquark is assumed to scatter as a single entity, 
p = (q,d). Note that CL is below 0.1%, so the quality of this fit is not acceptable. The bottom panel indicates, that the shape 
of the diffractive cross-section around the first diffractive minimun is not reproduced correctly by this model at LHC energies, 
and - as can also be seen on this Figure - this shortcoming cannot be fixed by leaving out a few data points around this dip 
from the optimalization procedure. 
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Proton=(qd), j4 qq and A fixed 



s = 23.5 GeV 



A qq =1.0 fixed 




FIG. 8. (Color online.) Visualisation of the obtained R q d, R q , Rd parameters for the case when the proton is assumed to 
scatter as a quark-diquark composit, p — (q, d) and A qq = 1, A = 0.5 fixed. The main observation is that the proton seems to 
be much larger at LHC energies than at ISR regime. This is mainly due to an increase in the R q d parameter, that characterizes 
the separation of the quark and the diquark. 
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[GeV] 


23.5 


30.7 


52.8 


62.5 


7000 


A 


0.5 


0.5 


0.5 


0.5 


0.5 


Aqq 


1.00 


1.00 


1.00 


1.00 


1.00 


R w ffml 


0.31 ± 0.01 


0.34 ± 0.01 


0.36 ± 0.01 


0.31 ± 0.01 


0.52 ± 0.01 


Ra ffml 


0.26 ± 0.01 


0.25 ± 0.01 


0.25 ± 0.01 


0.28 ± 0.01 


0.39 ± 0.01 


Rd [fm] 


0.71 ± 0.01 


0.71 ± 0.01 


0.71 ± 0.01 


0.75 ± 0.01 


0.83 ± 0.01 


X 2 /NDF 


65.0/44 


26.7/32 


57.2/33 


52.9/33 


4622.7/76 


CL [%] 


2.14 


73.05 


0.55 


1.56 


0.0 


^elastic [mb] 


6.2 ± 0.1 


5.1 ± 0.3 


5.0 ± 0.3 


7.3 ± 0.3 


17.9 ± 0.2 



TABLE I. The overall fit quality and resulting parameters of the fit at ISR energies including the LHC result at 7 TeV. The 
diquark is assumed to be a single entity. 

B. Model p= (q,(q,q)): diquark scatters as composite object 

In this subsection, similar MINUIT fit results are presented as in the previous subsection, the main modification 
is a change in the model assumption: now we assume that the diquark scatters as a composite object that contains 
two quarks. We present the final fit results to proton-proton elastic scattering data both at ISR [Ifl [26[ and at LHC 
energies Q- The results are illustrated on Fig. [9lfT3l The confidence levels, and model parameters with their errors 
are summarized in Table [TTJ 
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p+p -» p+p, diquark as a qq entity at \ s=23.5 GeV 
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FIG. 9. (Color online.) Same as Fig. [3] but the diquark is assumed to have a composite (q, q) substructure. 
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p+p p+p, diquark as a qq entity at \Js=30.7 GeV 
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FIG. 10. (Color online.) Same as FigH but for y/s = 30.7 GeV. 
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p+p -» p+p, diquark as a qq entity at \ s=52.8 GeV 
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FIG. 11. (Color online.) Same as FigH but for y/s = 52.8 GeV. 
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p+p -» p+p, diquark as a qq entity at \ s=62.5 GeV 
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FIG. 12. (Color online.) Same as FigEl but for = 62.5 GeV. 
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p+p -» p+p, diquark as a qq entity at \Js= 7000.0 GeV 




FIG. 13. (Color online.) Same as Fig[9] but for i/s = 7 TeV. Note that this fit is not acceptable, since its CL is below 0.1%. 



C. Total cross sections in the p = (q,(q,q)) model: the diquark is considered to be a composite object. 

The total inelastic cross sections for the quark-quark, quark-diquark and diquark-diquark subcollisions were analyzed 
on the basis of formula The detailed results are collected in Table Hill while the average of these ratios for the 
described ISR energies are 



20 



[GeV] 


23.5 


30.6 


52.8 


62.5 


7000 


A 


0.50 


0.50 


0.50 


0.50 


0.50 


Aqq 


1.00 


1.00 


1.00 


1.00 


1.00 


if w ffml 

±L qd L J 


0.21 ± 0.01 


0.23 ± 0.01 


0.25 ± 0.01 


0.23 ± 0.01 


0.66 ± 0.01 


Ra ffml 


0.28 ± 0.01 


0.28 ± 0.01 


0.28 ± 0.01 


0.29 ± 0.01 


0.33 ± 0.01 


Rd [fm] 


0.72 ± 0.01 


0.74 ± 0.01 


0.76 ± 0.01 


0.76 ± 0.01 


0.67 ± 0.01 


X 2 /NDF 


80.0/44 


32.1/32 


65.4/33 


68.6/33 


4852.7/76 


CL [%] 


0.07 


46.30 


0.07 


0.03 


0.0 


^elastic [mb] 


6.9 ± 0.1 


6.6 ± 0.1 


6.6 ± 0.1 


7.2 ± 0.1 


9.8 ± 0.1 



TABLE II. The overall fit quality and resulting parameters of the fit at the ISR energies including the LHC result at 7 TeV. 
The diquark is assumed to be a d — (q, q) composit entity. 



[GeV] 


23.5 


30.6 


52.9 


62.5 


7000 


a ~qd/ IJ qq 


1.92 ± 0.01 


1.93 ± 0.01 


1.93 ± 0.01 


1.93 ± 0.01 


1.88 ± 0.01 


°"dd/°"<jq 


3.64 ± 0.02 


3.66 ± 0.01 


3.67 ± 0.01 


3.65 ± 0.01 


3.43 ± 0.02 



TABLE III. The ratios of the total inelastic cross sections for the quark-quark, quark-diquark and diquark-diquark processes 
for the ISR and LHC energies using the composite diquark hypothesis. 



<J qq : <J qd :<r dd = l: (1.93 ± 0.03) : (3.65 ± 0.1), (18) 

which is close to the ideal 1:2:4 ratio, confirming the assumption of having two quarks inside the diquark, amended 
with some shadowing which is 4% and 9% respectively. At 7 TeV the ratios are different from ([18]) 

1 : (1.88 ± 0.01) : (3.43 ± 0.02), (19) 

which shows that shadowing is stronger, 6% and 14% percent respectively, and a significant decrease compared to the 
ideal ratio can be observed. 

The total elastic scattering cross sections were also determined, as given in Table [Til In a preliminary conference 
proceedings [j| , we have noted that that the errors on the total elastic cross sections might be decreased by fixing the 
A qq and A parameters to their nominal values of 1 and 0.5, respectively. Our final results presented here confirm this 
conjecture. 



[GeV] 
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30.6 


52.9 


62.5 


7000 


data [mb] 


38.94 ± 0.17 


40.14 ± 0.17 


42.67 ± 0.19 


43.32 ± 0.23 
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p=(q>d) 


vtotat [mb] 


38.5 ± 0.2 
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42.5 ± 0.2 


43.2 ± 0.3 


79.0 ± 1.1 


X 2 / NDF 


107.2/45 


62.2/33 


110.7/34 


56.1/34 


4667.67/77 


CL [%] 


0.00 


0.16 


0.00 


1.00 


0.00 


p=(q,(q,q)) 


vtotal [mb] 


38.6 ± 0.1 


39.7 ± 0.2 


41.6 ± 0.2 


42.6 ± 0.2 


74.2 ± 1.3 


X 2 / NDF 


86.8/45 


60.2/33 


154.2/34 


95.3/34 


5059.57/77 


CL [%] 


0.02 


0.26 


0.00 


0.00 


0.00 



TABLE IV. Measured and fitted values of the total cross sections, when the measured values of the total cross sections are 
included into the earlier fitting procedure, as another, additional data point. 
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Proton=(q(qq)), A qq and A fixed 
VT= 23.5 GeV 
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FIG. 14. (Color online.) Visualisation of the obtained R q d, R q , Rd parameters when the diquark is assumed to be a qq entity. 
The proton seems to be much larger at LHC energies than in the ISR regime. This is mainly due to an increase in the R q d 
parameter, that characterizes the separation of the quark and the diquark, which is further decomposed as (q,q) state in this 
class of models. 
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D. Model comparison 



The -y/s dependence of the parameters is summarized here. The confidence levels as a function of the center of mass 
energy for both of the models are given on Fig. Qj)l We left out the confidence level for TOTEM which is zero, all 
the other values are higher than or close to 0.1%. 

The center of mass energy dependence of the quark-diquark distance is presented on Fig. 1161 The effective quark 
radius can be seen on Fig. [T7] and the obtained diquark radius is on Fig. [TH1 

As can be seen on these figures, the fit quality is similar and the best fit parameters are rather different for both 
models at each colliding energies and the picture of the protons as realized by Figs. 151 and [TH are rather different. 
We started to wonder how it is possible, that the apparently rather different p=(q,d) and p=(q,(q,q)) models give so 
similar fit results, not only quantitatively but qualitatively too? Although the mathematical equations that describe 
the da/dt from the two models are rather complex and formally different from one another, we investigated some 
simple relationships among the model parameters. 

Numerically, we found one combination of the model parameters, an effective radius that is obtained from the 
quadratic sum of R q , Rd and R q d, that seems to be the same in both models and which is related in a simple way to 
the measured total cross-section. 




(20) 



G total = 27T-R f f . 



(21) 



The y/s dependence of the effective radius R c g is shown on Figure HH and the relationship between the effective radius 
and the measured total cross section is presented on Figure [501 
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The confidence level as a function of s 
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FIG. 15. (Color online.) The confidence level as a function of y/s. The value of CL at 7 TeV is practically zero and it is not 
shown. 



24 



The quark-diquark distance as a function of s 
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FIG. 16. (Color online.) Parameter R q d, representing the separation of the quark and the diquark, as a function of y/s. 
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The quark radius as a function of s 
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FIG. 17. (Color online.) Parameter R q (quark size) as a function of \/s. 
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The diquark radius as a function of s 
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FIG. 18. (Color online.) Parameter Rd, the diquark size, as a function of *Js. 
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The effective radius R eff as a function of s 
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FIG. 19. (Color online.) The model independent R c g radius parameter as a function of *Js. 
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FIG. 20. (Color online.) The effective i? c ff radius divided by the measured total cross section as a function of *Js. The result 
indicates a simple relationship between these parameters, namely, both for the case when a proton is modelled as p = (q,d) 
and for the case when p = (q, (q,q)), the quadratic sum of the characteristic radii is the same model independent quantity. 
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IV. CONCLUSION AND OUTLOOK 

A systematic study of fit quality as well as the fit parameters under similar circumstances have been performed for 
the Bialas - Bzdak model [l[ in a wide energy range from ISR to LHC energies using the same kinematic interval and 
the same method at each energy. The model gives a good description of the ISR data, which means that the CL is 
acceptable on the ISR energies if 3 data points at first diffractive minimum were left out from the fit. These studies 
confirm earlier results [l| concerning the increase of the total proton-proton cross section ( "size" of the proton) with 
increased colliding energies. 

An important shortcoming of the quark-diquark model of protons is that it ignores the real part of the elastic 
scattering amplitude. This leads to a singular behaviour at the diffractive minimum, which is apparently a more and 
more serious model limitation with increasing the energy of p+p collisions. Due to this reason, the considered model 
fails to describe in detail the structure of the first diffractive minimum at the LHC energies of 7 TeV. 

We found a combination of the model parameters, an effective radius that is obtained as the quadratic sum of quark, 
diquark radii and the separation between the quark and the center of mass of the diquark, that seems to be the same 
in both models and which is related in a simple, intuitive way to the measured total cross-section, as given by eqs. 
([2H)l and (|2"Tj) . As indicated on Figure [201 the precision of this "rule of thumb" formula is about 10 % which is quite 
amazing for us given that it is an extremely simple formula as compared to the full, exact and analytic expressions 
that are also obtained from both models. 

The evaluated ratios of the quark-quark, quark-diquark and diquark-diquark total inelastic cross-sections were 
found to deviate more and more from the ideal 1:2:4 ratio with increasing energies. In the ISR energy range the 
deviations from this ideal value were less than a 5<r effect, indicating lack of significant shadowing effects. However 
at the current LHC energy of y/s = 7 TeV, a significant decrease - as indicated by eq. (IT9)) - compared to these ideal 
ratios were found, which possibly may indicate an increased role of shadowing at CERN LHC energies. 

Finally let us note that the TOTEM Collaboration extended recently the measurement of the differential elastic 
p+p scattering cross-sections to low values of \t\ in ref. [27|, allowing one to extrapolate to the optical point at t — 
and to determine the total elastic and the total scattering cross-sections of p+p collisions at y/s = 7 TeV for the first 
time, but these data were yet not utilized in our analysis. 
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